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13.  AtSTRACT  (Minimum  20Q  wonts)  _ _ 

(1)  Selective  "blindness"  to  approaching  or  receding  motion  in  depth  exists  and  seems 
to  be  not  uncommon  in  normally-sighted  individuals.  Qf  16-subjects,  8  had  visual  field  » 
defects  for  either  approaching  or  receding  motion.  Of  21  subjects,  only  6  had  full  ' 
symmetrical  fields  for  oscillatory  motion  in  depth.  Visual  sensitivity  to'sideways  motion 
was  normal  in  stereomotion-blind  areas.  The  possible  relevance  to  aviation  is  pointed  out. 
.(2)  A  perfectly  camouflaged  bar  within  a  random  dot  pattern  was  rendered  visible  by 
moving  dots  within  the  bar  and  outside  the  bar  with  equal  and  opposite  velocities  (motion 
parallax).  The  feafVorientatioivcould  be  judged  with  equal  precision  (0.5  deg)  to  that  of  an 
uncamouflaged  dotted  bar  made  visible  by  brightness  contrast  providing  that  dot  speed  and 
contrast  were  high.  But  when  contrast  was  reduced,  discrimination  collapsed  for  the 
camouflaged  bar  (3)  Shape  discrimination  was  compared  for  motion-defined  and  contrast- 
defined  dotted  rectangles.  At  high  dot  speeds  and  contrasts  aspect  ratio  discrimination 
equal  for  the  two  kinds  of  rectangle  and,  at  2-3%,  corresponded  to  a  change  of  side  length 
of  only  24  arc  sec.  (4)  Orientation  discrimination  and  shape  discrimination  degrade  more 
rapidly  at  short  presentation  durations  for  a  motion-defined  than  for  a  contrast-defined 
target.  (5)  The  findings  in  (2)  -  (4)  above  suggest  that  helicopter  pilots  may  be  at  risk  of 
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making  visual  judgement  errors  in  nap  of  the  earth  flight  where  some  objects  and  ground 
features  are  seen  by  motion  alone  when  contrast  or  speed  is  low  or  when  inspection 
duration  is  brief.  (6)  We  have  developed  a  simple  portable  test  for  assessing  visual  ability 
to  see  and  discriminate  motion-defined  formrlt  is  aletter  readings  test  suitable  for-fieM  use  *' 
and  for  medical -use.  (7)  The  motion-defined  letter  test  was  used  on  25  patients  with 
multiple  sclerosis  and  50  controls;  34/50  eyes  of  patients  were  abnormal  even  though  visual 
acuity  was  normal.  (8)  Nonlinear  systems  analysis:  We  have  developed  a  new 
mathematical  approach  to  testing  multi-neuron  models  in  which  individuals  neurons  are 
modelled  as  rectifiers.  (9)  We  have  developed  a  nondestructive  zoom-FFT  technique  that 
allows  spectra  of  EEO  and  other  time  series  to  be  computed  with  the  theoretical  resolution 
allowed  by  the  Heisenberg-Gabor  relation,  e.g.  50,000  lines  DC- 100  Hz  at  a  resolution  of 

O. 002  Hz  from  a  500-sec  recording.\(10)  By  using  a  2-sinewave  nonlinear  analysis 
approach  in  recording  human  evoked  potentials  we  have  found  that  both  vertically- tuned 
and  horizontally-tuned  responses  have  a  bandwidth  of  about  12  deg,  and  that  there  is  a 
strong  nonlinear  interaction  between  horizontal  and  vertical.  (11)  Single-cell  responses: 

We  have  built  electronics  to  allow  single  cell  recordings  to  be  analyzed,  and  carried  out 
pilot  experiments.  (12)  A  magnetically-shielded  room,  7-channel  magnetometer  and 
sensory  position  indicator  (SPI)  have  been  installed  at  York  University  and  several  faults 
rectified.  We  have  investigated  nonlinearities  in  the  response  to  two  visual  stimuli  of 
different  frequencies  and  two  auditory  stimuli  of  different  frequencies,  analyzing  the 
responses  by  zoom-FFT.  (13)  A  book  "Human  Brain  Electrophysiology”  written  by  the 

P. I.  was  published  in  1989.  (14)  Two  books  edited  by  the  P.I.,  one  on  "Binocular 
Vision"  and  one  on  "Spatial  Vision  are  in  press  for  1990  publication. 
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2a.  Objectives:  Psychophysical 

(1)  Further  define  the  roles  of  the  channeling  hypothesis  in:  (a)  identifying  specific  visual 
processes;  (b)  understanding  visual  performance;  (c)  specifying  visual  parameters  likely  to  be 
important  in  eye-hand  coordination,  especially  in  aviation  and  flight  simulator  visual  displays. 

(2)  Camouflage  and  the  visual  processing  of  objects  defined  by  motion  alone.  For  camouflaged 
objects  that  are  invisible  except  when  there  is  motion  parallax  between  the  object  and  background, 
measure  spatial  discriminations,  and  in  particular  the  hyperacuities,  orientation  discrimination, 
spatial  frequency  discrimination,  and  line  interval  discrimination.  Compare  these  data  with  the 
corresponding  hyperacuities  for  objects  defined  by  luminance  contrast,  and  find  whether  both  sets 
of  data  can  be  explained  by  an  opponent  or  line-element  model  of  spatial  form  discrimination 
proposed  previously/1'5) 

2a.  Objectives:  Neuromagnetism  and  electrophysiology 

(1)  Link  the  channeling  modes  of  human  psychophysics  with  the  activation  of  different  sensory 
projection  areas  in  human  cortex. 

(2)  Identify  evoked  activity  in  different  visual,  auditory  or  somatosensory  projections  in  human 
cortex  and  elucidate  the  differences  between  the  type  of  processing  occurring  in  the  different  areas. 
Link  these  data  with  the  known  functional  neuroanatomy  of  macaque  monkey  brain  and  with 
human  psychophysics. 

(3)  Elucidate  the  temporal  sequence  of  activation  of  different  cortical  areas  evoked  by  different 
kinds  of  complex  visual  auditory  and  somatosensory  stimuli.  These  data  will  complement  scan¬ 
ning  data  (e.g.  regional  cerebral  blood  flow,  PET)  that  lack  the  temporal  resolution  offered  by 
neuromagnetic  recording. 
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(4)  Elucidate  relationships  between  simultaneous  activities  of  different  cortical  areas  within  a 
single  modality  (visual,  auditory  or  somatosensory). 

(5)  Identify  the  cortical  sites  of  interactions  between  responses  to  stimuli  of  different  modalities, 
and  compare  these  sites  with  the  known  poly-sensory  cortical  areas  in  nonhuman  primates. 

(6)  By  combining  neuromagnetic  and  evoked  potential  recording,  exploit  their  complementary 
natures  to  improve  the  localization  of  generator  sites. 

(7)  Locate  the  brain  sites  of  abnormalities  in  patients  with  known  specific  sensory  defects 
including  selective  orientation-tuned  visual  loss  for  intermediate  spatial  frequencies,  stereomotion 
"blindness",  specific  defects  of  shape  recognition,  selective  deafness  to  frequency  changes. 

2b.  Status  of  the  Research  Effort:  Psychophysical 

(1)  Specific  "blindness"  to  oscillatory  and  unidirectional  motion  in  depth. 

As  an  object  moves  towards  the  head  its  two  retinal  images  move  in  opposite  directions. 
This  binocular  cue  alone  can  generate  a  strong  impression  of  motion  in  depth  (stereomotion).  We 
have  previously  published  visual  fields  for  oscillatory  motion  in  depth  and  found  that  normally- 
sighted  subjects  have  areas  of  specific  blindness  to  stereomotion/6-7)  Of  the  six  subjects  reported, 
five  showed  stereomotion  field  defects.  We  have  now  extended  the  data  base  to  a  further  21 
normal  subjects,  and  confirm  that  stereomotion  field  defects  are  common.  Only  6/21  subjects  had 
full  symmetrical  fields. 

We  now  report  the  existence  of  selective  blindness  to  unidirectional  motion  in  depth.  Of  16 
subjects  whose  visual  fields  were  tested  for  approaching  and  for  receding  motion  in  depth,  only  5 
had  similar  fields  for  approaching  and  receding  motion. 
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Table  1  summarizes  the  data.  Figure  1  illustrates  stereomotion  fields  that  were  full  and 
symmetrical.  Figure  2  shows  fields  for  a  subject  with  field  defects  and  areas  that  were  "blind"  to 
motion  in  one  direction. 

Because  sensitivity  to  monocularly-viewed  motion  showed  no  abnormalities  corresponding 
to  the  binocular  stereomotion  "blind  spots"  we  conclude  that  the  stereomotion  field  defects  were 
chiefly  due  to  the  cortical  processing  of  motion.  We  also  conclude  that  unidirectional  motion 
defects  are  caused  by  a  loss  of  sensitivity  to  unidirectional  motion  in  depth  rather  than  to  abnormal 
interactions  between  mechanisms  for  approaching  and  receding  motion.  These  findings  provide 
further  evidence  that  the  human  visual  pathway  contains  different  binocular  mechanisms  for 
position  in  depth  and  for  motion  in  depth,  and  that  stereomotion  blindness  is  due  to  a  selective  loss 
of  the  motion  mechanism. 

These  findings  raise  the  possibility  that  stereomotion  "blind  spots"  are  not  uncommon  in 
pilots,  and  that  the  trajectory  of  an  oncoming  aircraft  might  be  misjudged  if  it  passed  through  a 
stereomotion  "blind  spot". 

A  report  on  the  results  to  date  has  been  published  in  Vision  Research.®) 

(2)  Orientation  discrimination  for  camouflaged  objects  defined  by  motion  alone  and  for  objects 

defined  by  luminance  contrast 

A  pseudo-random  pattern  of  bright  dots  subtending  2.2  x  2.2  deg  was  generated  by 
hardware  of  our  own  design.  Frame  rate  was  200  Hz.  Dots  subtended  2.0  min  arc,  mean 
separation  was  about  6  min  arc  and  there  were  approximately  1000  dots.  The  dots  were  optically 
superimposed  on  a  circular  uniformly-illuminated  area  of  diameter  3.7  deg.  A  1.5  x  0.22  deg  bar- 
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Table  1.  Summary  of  results  for  21  normally-sighted  subjects.  Key: 
D-  different  fields.  S  -  similar  fields.  U  -  unclassified.  U(PR)  - 
unclassified  with  poor  reproducibility. 
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Figure  1.  Visual  fields  for  unidirectional  depth  perception . 
A  subject  with  similar  large  fields  for  approaching  and 
receding  motion. 
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Figure  2.  Large  visual  field  defect  for  unidirectional  depth 
perception.  A,  B;  near  disparities,  approaching  motion. 

C,  D:  near  disparities,  receding  motion. 
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shaped  area  within  the  dot  pattern  was  rendered  visible  by  moving  dots  inside  this  area  leftwards 
and  outside  this  area  rightwards  at  constant  velocity.  When  the  dots  were  stationary  the  bar  was 
perfectly  camouflaged.  Dot  contrast  was  varied  by  neutral  density  filters.  Orientation  discrimina¬ 
tion  was  measured  by  temporal  two-alternative  forced  choice.  The  dot  pattern  was  presented  for 
1.0  sec,  and  contained  a  motion-defined  vertical  bar.  Then  there  was  an  interval  of  0.5  sec 
followed  by  a  second  presentation  of  1.0  sec  with  the  bar  inclined  at  some  angle  9.  There  were  10 
possible  values  of  0.  Bar  location  was  randomly  jittered  and  a  fresh  random  dot  pattern  was 
generated  for  each  presentation.  The  subject  pressed  one  of  two  buttons  depending  on  whether  9 
was  clockwise  or  anticlockwise  of  vertical.  Orientation  discrimination  threshold  was  calculated  by 
Probit  analysis. 

In  separate  experiments  orientation  discrimination  was  measured  for  a  non-camouflaged  bar 
that  was  created  by  omitting  the  dots  in  the  area  surrounding  the  bar.  This  target  is  illustrated  in 


Figure  3.  A  -  random  dot  pattern  containing  a  perfectly  camouflaged  bar.  B  -the  bar  was 
revealed  by  moving  dots  within  the  bar  and  outside  the  bar  in  opposite  directions. 
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The  rationale  of  this  experiment  is  that,  for  the  camouflaged  bar,  figure  ground  segregation 
was  achieved  by  motion  alone  but  the  non-camouflaged  bar  was  rendered  visible  by  luminance 
contrast  Dot  density  and  velocity  within  the  bar  were  identical  in  the  two  cases. 

Figures  4  and  5  show  that,  for  high  dot  velocities  and  contrasts,  orientation  discrimination 
is  similar  for  motion-defined  and  contrast-defined  bars.  Furthermore,  at  about  0.4  deg,  discrimi¬ 
nation  compares  favourably  with  the  most  acute  values  reported  in  the  literature  for  conventional 
bright  solid  bars  or  lines.  This  finding  may  relate  to  our  previous  finding  that  vernier  acuity  for  a 


camouflaged  dotted  bar  can  be  as  high  as  for  a  non-camouflaged  dotted  bar  (see  Final  Report  dated 


1987/09/14  and  Reference  #9). 
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Figure  4.  Orientation  discrimina¬ 
tion  versus  dot  speed  for  a  dotted  bar 
defined  by  relative  motion  (open 
symbols)  and  for  the  same  bar 
defined  by  luminance  contrast  (filled 
symbols).  Bar  detection  thresholds 
are  arrowed 
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3.5  r 


CD 


DOT  CONTRAST 

Figure  5.  Orientation  discrimination  versus  dot  contrast.  Other  details  as  in  Figure  4. 
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But  Figure  5  also  shows  that,  as  contrast  is  reduced,  discrimination  collapses  earlier  for  the 
motion-defined  bar  than  for  the  contrast-defined  bar.  In  particular,  there  is  a  contrast  range  of 
about  4: 1  over  which  discrimination  has  collapsed  for  the  motion-defined  bar  but  is  still  good  for 
the  contrast-defined  bar. 


ABC 


Figure  6.  A-C  are  three  snapshots  of  the  dotted  bar  taken  during  a 
1.0  sec  presentation.  The  dots  surrounding  the  bar  were  switched  off. 
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Turning  back  to  Figure  4,  we  now  consider  the  effect  of  dot  velocity  on  discrimination  for 
the  camouflaged  motion-visible  bar.  It  is,  in  principle,  possible  that  dot  motion  might  improve 
discrimination  by  reducing  errors  due  to  spatial  sampling  by  dots.  Figure  6  illustrates  this  point. 
Because  of  the  coarse  spatial  sampling  provided  by  the  sparse  dots,  the  orientation  of  the  bar's 
edge  is  poorly  defined  in  each  photograph.  But  in  principle,  orientation  could  be  more  precisely 
defined  by  taking  all  three  "snapshots"  into  account.  However,  Figure  4  shows  that  this  effect  did 
not  occur  for  the  contrast  defined  bar  (filled  symbols).  We  can  therefore  assume  that  the  effect  of 
velocity  on  discrimination  for  the  camouflaged  bar  (open  symbols)  was  due  to  velocity  sensitivity 
of  moiton-sensitive  mechanisms  rather  than  to  sampling  errors. 

A  preliminary  report  of  this  study  has  been  published  in  Vision  Research/10). 

(3)  Shape  discrimination  for  camouflaged  objects  defined  by  motion  alone  and  for  objects 
defined  by  luminance  contrast 

We  have  used  a  similar  technique  to  that  described  in  #2  above  to  generate  a  camouflaged 
rectangular  shape  that  is  visible  by  motion  alone.  The  rectangular  area  was  rendered  visible  by 
moving  the  dots  inside  the  rectangle  obliquely  downwards  to  the  left  at  a  constant  speed  that  was 
equal  and  opposite  to  the  speed  of  the  dots  in  the  remainder  of  the  pattern.  The  reason  for 
choosing  oblique  motion  was  to  ensure  that  the  rectangle's  vertical  and  horizontal  edges  were  all 
defined  by  the  same  mix  of  shearing  and  compressive  motion.  Although  the  rectangle  was  visible 
to  the  eye,  a  photographic  time  exposure  did  not  reveal  it  (Figure  7B).  Contrast-defined  rectangles 
were  created  by  switching  off  all  dots  outside  the  rectangle  (Figure  7C).  The  ratio  between  the 
lengths  of  vertical  and  horizontal  sides  had  10  possible  values,  and  these  were  presented  randomly. 
The  subject's  task  was  to  press  one  of  two  buttons  depending  on  whether  the  longer  sides  were 
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vertical  or  horizontal.  To  ensure  that  both  dimensions  must  be  compared,  different  areas  of 


rectangle  were  interleaved  randomly  as  illustrated  in  Figure  8.  To  ensure  that  the  distance  of  any 
edge  from  the  boundary  of  the  display  provides  no  cue  to  shape,  the  rectangle's  location  was 
jittered  randomly.  Shape  discrimination  threshold  was  measured  by  two-alternative  forced  choice 


and  Probit  analysis. 

3.7° 


Figure  7.  A-  Photograph  of  a  typical 
random  dot  pattern.  The  test  rectangle 
is  perfectly  camouflaged.  The  outer 
circle  marks  the  edge  of  the 
superimposed  adapting  field.  B -  Dots 
within  the  rectangle  and  outside  the 
rectangle  were  moved  at  equal  and 
opposite  speeds.  Although  the 
rectangle  was  evident  to  the  eye,  it  was 
not  evident  to  the  camera.  C-  All  dots 
outside  the  rectangle  were  switched  off 
to  create  oan  uncamouflaged  rectangle. 
Aspect  ratio  was  defined  as  (alb). 
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DISSOCIATION  OF  ASPECT  RATIO  (a/b) 
FROM  LINEAR  DIMENSIONS 

===  -  ■>% 


1.2 

AREA 

1.0 

AREA 


1/1.2 

AREA 


Figure  8.  Dissocation  of  aspect  ratio  (alb)  from  area  and  from  any  given  linear  dimension 
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Figure  9.  Ordinate  plot  shape  discrimination  thresholds  A  (alb)  expressed  as  a  percentage, 
where  A  (alb)  is  the  smallest  disriminable  change  in  the  ratio  A  (alb).  Abscissae  plot  dot 
speed  (relative  speed  is  twice  dot  speed).  Open  circles  are  for  motion-defined  rectangles  and 
filled  circles  are  for  uncamouflaged  rectangles.  Dot  contrast  in  B  and  D  was  0.6  log  units  less 
than  in  A  and  C.  The  leftmost  filled  circles  in  each  panel  is  for  zero  speed.  A,  B  -  subject  1. 
C,  D  -  subject  2.  Viewing  was  monocular. 
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We  have  measured  shape  discrimination  as  a  function  of  dot  speed  and  dot  contrast  for 
camouflaged  dotted  rectangles  and  for  uncamouflaged  dotted  rectangles.  Figure  9  shows  that, 
when  dot  speed  and  contrast  were  both  high,  aspect  ratio  discrimination  threshold  was  as  acute  for 
a  motion-defined  rectangle  as  for  a  contrast-defined  rectangle  and,  at  2-3%,  corresponded  to  a 
change  of  side  length  of  about  24  arc  sec.  This  is  a  remarkable  visual  performance,  given  that  the 
mean  dot  separation  was  about  360  arc  sec.  Discrimination  collapsed  at  low  dot  speeds  and  could 
not  be  measured  at  speeds  less  than  about  0. 1-0.2  deg/sec  for  motion-defined  rectangles,  but  was 
almost  unaffected  by  dot  speed  for  contrast-defined  rectangles.  To  explain  how  the  visual  system 
can  dissociate  changes  of  shape  from  changes  of  area  we  invoke  a  previous  suggestion^1)  that  the 
visual  system  contains  a  mechanism  that  compares  the  separation  of  pairs  of  contours  along 
perpendicular  azimuths.  To  explain  the  coincidence  of  shape  discrimination  thresholds  for  motion- 
defined  and  contrast-defined  rectangles  we  suggest,  that  in  early  visual  development  the  underlying 
neural  mechanisms  are  driven  by  the  same  environmental  and  behavioural  pressures  towards  a 
common  endpoint.  In  other  words,  "spatial  vision  is  spatial  vision"  and  the  performance  required 
for  eye-hand  coordination  that  is  dictated  by  the  outside  environment  independently  of  how  the  eye 
segregates  figure  from  ground.  This  idea  could  also  explain  why  orientation  discrimination  is  the 
same  for  motion-defined  and  contrast-defined  bars. 

A  report  on  this  work  has  been  submitted  for  publication^2)- 
(4)  Dynamics  of  orientation  discrimination  and  shape  discrimination  for  objects  defined  by 
motion  along  and  for  objects  defined  by  luminance  contrast 
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We  have  studied  the  dynamics  of  orientation  discrimination  and  shape  discrimination  in  the 
following  two  ways:  (A)  By  measuring  discrimination  threshold  as  a  function  of  presentation 
duration,  dot  lifetime  being  equal  to  presentation  duration;  (B)  by  measuring  discrimination 
threshold  as  a  function  of  dot  lifetime  for  a  fixed  presentation  duration.  These  two  studies  are  in 
progress.  Results  to  date  are  as  follows. 

Figure  10  shows  that,  for  presentation  durations  shorter  than  about  0.15  sec,  shape 
discrimination  rapidly  became  poorer  for  MD  than  for  CD  rectangles.  This  effect  could  not  be 
attributed  entirely  to  an  elevation  of  simple  motion  detection  threshold  at  low  durations  because  the 
effect  of  duration  on  motion  threshold  (Figure  10B)  was  quite  different  to  the  effect  of  duration  on 
shape  discrimination  (Figure  10A).  In  terms  of  a  Reichardt-type  model  of  figure-ground 
discrimination  we  might  regard  Figure  10B  as  reflecting  the  properties  of  peripheral  motion 
detectors,  and  Figure  10B  as  reflecting  the  properties  of  figure-ground  extraction  and  more  central 
processing.  In  these  terms,  the  difference  curve  (Figure  IOC)  reflects  the  properties  of  figure- 
ground  extraction  and  more  central  processing. 

The  effect  of  presentation  duration  on  orientation  discrimination  is  similar  to  the  effect  on 
shape  discrimination  shown  in  Figure  10A.  The  effect  of  dot  lifetime  on  orientation  discrimination 
is  somewhat  different.  Discrimination  improves  with  dot  lifetime  up  to  about  100  msec,  beyond 
which  it  has  little  effect. 
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Figure  10.  A-  The  effect  of 
presentation  duration  on  shape 
discrimination  for  camouflaged 
motion-defined  rectangles.  Just- 
discriminable  differences  in  aspect 
ratio  (ordinate)  are  plotted  versus 
duration  (abscissa)  for  5  different  dot 
speeds.  Open  square,  plusses,  filled 
stars,  open  circles,  and  filled  circles 
respectively  plot  dot  speeds  of  0.18, 
0.36,  0.68,  1.36  and  2.72  deglsec. 
Aspect  ratio  is  defined  in  Fig.  ID. 
The  velocity  difference  between  dots 
in  the  rectangle  and  surround  equals 
twice  the  dot  speed.  B-  The  effect  of 
presentation  duration  (abscissa)  on 
velocity  direction  discrimination 
threshold  defined  as  the  lowest  dot 
speed  for  which  leftward  and 
rightward  motion  could  be 
discriminated  with  75%  accuracy. 
Both  axes  are  logarithmic  in  A  and 
B.  C-  Difference  between  the  curves 
in  A  and  B.  Subject  #/. 
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Reports  on  the  work  completed  to  date  have  been  submitted  for  publication^3-14). 

(5)  A  simple  device  for  testing  visual  sensitivity  to  motion-defined  objects 

From  the  viewpoint  of  aviation,  the  significance  of  the  findings  reported  in 
Section  2-5  above  is  that  they  suggest  that  in  nap  of  the  earth  helicopter  flight, 
where  some  ground  features  are  visible  by  motion  alone  while  others  are  visible 
by  contrast,  a  pilot's  visual  judgements  might  fail  for  motion-visible  objects  but 
not  for  contrast-visible  objects  even  though  the  motion-visible  objects  are  still 
clearly  visible.  From  the  viewpoint  of  medical  research,  the  significance  is  the 
suggestion  that  spatial  vision  for  motion-defined  form  is  mediated  by  different 
neurons  than  spatial  vision  for  contrast-defined  form  so  that  there  may  be  forms 
of  visual  damage  or  visual  abnormality  that  do  not  affect  visual  acuity  or  contrast 
sensitivity,  but  degrade  spatial  vision  for  motion-defined  objects. 

It  is  not  practical  to  follow  up  these  implications  with  the  equipment  used  to  obtain  the  data 
reported  above,  because  that  is  "one-off"  equipment  of  our  own  design  and  is,  complex  and 
nonportable.  Therefore  we  have  developed  a  simple  device  for  testing  visual  sensitivity  to  motion- 
defined  form  with  the  intent  of  providing  a  means  of  applying  this  basic  research  in  field  studies  of 
pilots  and  in  hospital  medical  research. 
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DOT  SPEED  (deg /sec) 


Figure  12  Ordinates  plot  visual  acuity  expressed  as  a  Snellen  fraction  con  a  log 
scale.  Abscissae  plot  dot  speed  on  a  linear  scale.  Open  symbols  indicate  that 
letters  were  rendered  visible  by  impressing  equal  and  opposite  speeds  on  the 
dots  inside  and  outside  the  letter.  Filled  symbols  Indicate  that  letters  were 
uncamouflaged  as  in  Fig.  1C.  Dot  contrast  was  near  1.0  and  dot  density  was 
20%.  Presentation  duration  was  4  sec  viewing  distance  was  9  m.  Panels  A  -  C 
respectively  show  data  for  subjects  1,  3  and  5. 
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Figure  1 1  illustrates  the  principle.  Figure  1 1 A  is  a  photograph  of  a  random  dot  pattern 
generated  on  the  screen  of  standard  IBM  PC.  A  letter  (Z)  at  the  centre  of  the  screen  is  perfectly 
camouflaged.  In  FigurellB,  dots  within  the  Z  move  rightwards  and  dots  outside  the  Z  move 
leftwards.  The  Z  is  visible  to  the  normal  eye,  but  is  not  visible  to  the  camera.  In  Figure  1 1C,  all 
dots  outside  the  Z  have  been  switched  off  so  that  the  Z  is  defined  by  luminance  contrast.  The 
subject  is  required  to  read  the  letters  presented.  There  are  10  letters  for  each  of  9  letter  sizes. 
Visual  acuity  for  motion-defined  letters  can  be  measured  at  some  given  dot  speed  by  plotting 
percent  correct  reading  score  versus  letter  size  on  probability  paper,  and  reading  off  the  letter  size 
for  a  75%  correct  score.  Alternatively,  speed  threshold  for  a  given  letter  size  by  plotting  percent 
correct  reading  score  versus  dot  speed  on  probability  paper,  and  reading  off  the  dot  speed  for  a 
75%  correct  score. 

In  order  to  thoroughly  "de-bug"  this  test,  an  extensive  parametric  study  on  normally- 
sighted  subjects  was  undertaken,  and  this  has  been  published^15).  Figure  12  (open  circles) 
shows  that  visual  acuity  for  motion-defined  letters  falls  off  slowly  as  dot  speed  is  reduced 
below  0.3  deg/sec  until  speed  falls  below  a  critical  value  below  which  acuity  cannot  be 
measured.  On  the  other  hand,  visual  acuity  for  a  contrast-defined  letter  (Figure  11 C)  is 
essentially  independent  of  dot  speed.  Open  and  closed  circles  in  Figure  13  show, 
respectively,  the  effects  of  dot  contrast  on  visual  acuity  for  motion-defined  and  contrast- 
defined  dotted  letters.  The  effect  of  dot  density  is  shown  in  Figure  14.  Acuity  for  motion- 
defined  letters  (open  circles)  was  comparatively  unaffected  by  dot  density  from  50%  to 
about  0.5%,  (a  1000:  1  range)  below  which  it  abruptly  collapsed.  On  the  other  hand,  acuity 
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Figure  13.  Ordinates  plot  visual  acuity  expressed  as  a  Snellen  fraction  for 
dotted  letters  on  a  log  scale.  Abscissae  plot  dot  contrasts  on  a  log  scale.  Open 
symbols  indicate  that  letters  were  rendered  visible  by  impressing  equal  and 
opposite  speeds  of  0.47  deg/sec  on  dots  inside  and  outside  the  letter.  Filled 
symbols  indicate  that  letters  were  uncamouflaged  as  in  Fig.  1C  with  dot  speed 
of  0.47  deglsec.  Dot  density  was  20%.  Viewing  distance  was  6m.  Presentation 
duration  was  4  sec.  Panels  A-D  respectively  show  data  for  subjects  1,  3,  2  and 
6. 
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DOT  DENSITY  (percent) 


Figure  14  Ordinates  plot  visual  acuity  expressed  as  a  Snellen  fraction  for 
dotted  letters  on  a  log  scale.  Abscissae  plot  the  density  of  bright  dots  on  a 
log  scale.  Open  symbols,  letters  rendered  visible  by  impressing  equal 
and  opposite  speeds  of  0.31  deg/sec  on  dots  inside  and  outside  the  letters. 
Filled  symbols ,  uncamouflaged  letters  as  illustrated  in  Fig.  1C  with  dot 
speed  of  0.31  deg/sec.  Presentation  duration  was  4  sec.  Dot  contrast  was 
/.0.  Viewing  distance  was  9m,  Panels  A-Z)  respectively  show  data  for 
subjects  1,  3,  5  and  4. 
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for  contrast-defined  letters  (filled  circles)  progressively  fell  as  dot  density  was  reduced  from  50%, 
and  below  about  0.5%  approximated  acuity  for  motion-defined  letters. 

(6)  Investigation  of  visual  loss  in  patients  with  multiple  sclerosis  using  the  motion-defined, 
letter  test  described  in  Section  6 

We  have  investigated  50  eyes  of  25  patients  with  multiple  sclerosis,  and  50  eyes  of  50 
control  subjects  using  the  following  tests:  (A)  Speed  threshold  for  75%  correct  reading  score  on 
large  (6/60)  motion-defined  letters;  (B)  Contrast  threshold  for  75%  correct  reading  score  for  6/60 
contrast-defined  dotted  letter;  (C)  Visual  acuity  for  75%  correct  reading  score  for  motion-defined 
letters  with  a  fixed,  high,  dot  speed;  (D)  Visual  acuity  for  solid,  contrast-defined  letters  on  the 
Regan  low-contrast  letter  charts  (96%,  11%  and  4%  letter  contrasts);  (E)  Isolated-letter  visual 
acuity  using  the  same  letters  as  on  the  Regan  letter  charts. 

We  found  that  speed  threshold  for  motion-defined  letters  could  detect  visual  abnormality 
that  was  not  detected  by  any  of  the  other  tests  used.  Of  the  50  eyes  of  patients,  34  (i.e.  68%)  were 
abnormal  at  the  2.5  SD  level,  and  of  these  34  eyes,  5  were  effectively  motion  blind  in  the  sense 
that  they  could  not  read  large  letters  even  at  our  highest  speed  difference  of  0.9  deg/sec.  (Note  that 
this  is  a  new  kind  of  motion  blindness,  quite  different  from  the  motion-in-depth  blindness  reported 
above).  We  conclude  that  our  motion-defined  letter  test  can  detect  pathology  that  is  not  picked  up 
by  testing  contrast  sensitivity  or  low-contrast  acuity.  We  suggest  that  this  new  test  can  detect 
dysfunction  in  the  human  equivalent  of  a  pathway  in  monkey  brain  that  originates  in  large  retinal 
ganglion  cells,  passes  through  the  magnocellular  layers  of  the  lateral  geniculate  body,  includes 
cortical  area  MT,  and  is  involved  in  processing  motion. 
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A  manuscript  is  about  to  be  submitted  for  publication^9)  and  two  abstracts  are  submitted 

(20,21). 

(7)  Investigation  of  visual  loss  in  patients  with  ocular  hypertension  and  early  glaucoma  using 
the  motion-defined  letter  test  described  in  Section  6. 

This  study  is  in  progress.  Twenty  patients  have  been  tested. 

(8)  Investigation  of  visual  loss  in  patients  with  Parkinson's  disease. 

Toronto  Western  Hospital  is  one  of  the  sites  of  the  US/Canadian  multi-centre  trial  of  the 
new  drug  depranyl  and  Professor  A  Lang,  MD  is  Director  of  this  site.  This  drug  is  thought  to 
slow  the  progress  of  brain  (nigral)  damage  caused  by  the  disease,  though  it  does  not  control  the 
symptoms.  (While  conventional  therapy  controls  the  symptoms,  the  disease  continues  to 
progress).  We  and  others  have  previously  reported  evidence  suggesting  that  the  motion  pathway 
can  be  selectively  affected  in  Parkinson's  disease06-18).  jn  a  joint  study  with  Dr.  Lang  we  are 
asking  whether  our  motion-defined  letter  test  can  provide  an  index  of  the  functional  integrity  of 
dopaminergic  neurons  in  the  visual  system,  and  if  so  whether  this  also  provides  an  index  of  the 
patient's  general  clinical  condition. 

By  comparing  eye  movement  records  with  MR  images,  Dr.  J.A.  Sharp,  Head  of 
Neurology  at  the  University  of  Toronto  has  discovered  that  eye  movements  in  response  to  a 
moving  target  are  defective  in  humans  with  focal  brain  lesions  in  an  area  of  cortex  suspected  to  be 
equivalent  to  MT  and  MST  in  macaque.  It  is  known  that  damage  to  MT  in  macaque  causes  eye 
movement  defects  and  also  behavioural  defects  in  motion  processing.  We  are  currently  carrying 
our  a  joint  study  with  Dr.  Lang  to  compare  pursuit  eye  movements  with  reading  scores  for  motion- 
defined  targets  in  patients  with  Parkinson's  disease  and  patients  with  focal  brain  lesions. 
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2b.  Status  of  the  research  effort:  neuromagnetism  and  electro¬ 
physiology 

C 9)Theoretical  and  technical  work  on  the  two-input  technique  for  characterizing  nonlinear 
processing  in  sensory  pathways 

Neurons  that  respond  asymmetrically  -  e.g.  to  leftwards  versus  rightwards  motion,  to 
increase  versus  decrease  of  spatial  contrast,  or  to  rise  versus  fall  of  auditory  tone  frequency  -  can 
be  described  as  rectifiers.  In  addition  to  asymmetric  response,  many  neurons  perform  functionally 
-  important  nonlinear  processing  such  as  ratio-ing/22-24)  multiplication/25)  or  logarithmic 
compression. 

We  have  developed  a  theoretical  basis  and  a  practical  technique  for  investigating  nonlinear 
processing  in  sensory  pathways.  The  basic  procedure  can  be  traced  back  at  least  to  Bennet's  1933 
paper(26)  on  radio  communication.  In  general  terms,  Bennet's  basic  idea  was  to  stimulate  the 
nonlinearity  being  studied  with  two  simultaneous  inputs,  one  of  temporal  frequency  Fj  Hz  and  the 
other  of  F2  Hz.  Any  other  frequency  terms  must  be  due  to  nonlinear  processing. 

Benne/26)  discussed  the  case  of  simple  linear  rectifier,  and  showed  theoretically  that  the 
output  included  many  terms  of  frequency  (nFi  ±  1T1F2),  where  n  and  m  are  integral  or  zero. 

Bennet  considered  the  case  that  the  amplitude  of  the  Fj  Hz  input  is  held  constant  while  the 
amplitude  of  the  F2  Hz  input  is  progressively  increased  from  zero,  and  developed  a  theoretical 
method  for  calculating  how  the  amplitudes  of  the  several  discrete  frequency  terms  vary  with  the  Fj 
Hz  input  amplitude. 

Bennet's  theoretical  work  was  further  developed  by  Rice  but  was  not  extended  previously 
to  rectifiers  of  any  given  characteristic  nor  to  cascades  of  rectifiers. 
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We  have  made  the  following  further  steps.  We  have  developed  a  theoretical  treatment  of 
the  following  cases:  (a)  single  compressive  rectifier,  y  =  x1^;  (b)  single  accelerating  rectifier 
y  =  xn;  (c)  cascaded  sequence  of  rectifiers,  e.g.  multiple  compressive  third  root  rectifiers  in  series, 
and  mixed  cascaded  rectifiers,  e.g.  compressive  third  root  followed  by  accelerating  square  law; 
(d)  two  parallel  rectifiers  (compressive  or  accelerating)  converging  onto  a  third  (compressive  or 
accelerating);  (e)  a  single  rectifier  whose  characteristic  matches  the  physiological  contrast 
sensitivity  characteristics,  i.e.  a  threshold  -  initial  acceleration  -  subsequent  compression/27) 

A  sequence  of  cascaded  rectifiers  (c  above)  is  intended  to  model  a  sequence  of  rectifier-like 
neurons  as,  for  example,  the  photoreceptor-bipolar-ganglion  cell-LGN  cell-cortical  cell  sequence. 
Case  (d)  above  is  intended  to  model  the  dichoptic  visual  situation  (i.e.  signals  leaving  nonlinear 
processors  in  left  and  right  eyes  converging  onto  binocular  cortical  neurons)  or  the  dichotic 
situation  (i.e.  signals  leaving  nonlinear  processors  in  left  and  right  ears  converging  onto  binaurally- 
driven  cells). 

We  went  on  to  compute  the  amplitudes  of  several  (up  to  20)  of  the  discrete  nonlinear 
frequency  components  as  a  function  of  the  amplitude  of  the  F2  Hz  input/27) 

In  brief,  this  theoretical  work  suggests  that  the  resulting  family  of  curves  comprises  a 
"fingerprint"  of  the  type  of  nonlinearity.  Because  so  many  different  frequency  components  are 
computed,  just  as  with  a  human  "fingerprint,"  there  is  high  specificity,  allowing  different  kinds  of 
nonlinearity  to  be  recognized. 

The  following  is  an  outline  of  this  mathematical  work.  A  full  treatment  of  the  work  to  date 
has  been  published  in  the  Journal  of  Theoretical  Biology Z28) 
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A  METHOD  FOR  DERIVING  THE  RESPONSE  OF  ASYMMETRIC  NONLIN¬ 
EARITIES  TO  A  SUM  OF  TWO  SINEWAVES 

We  first  consider  the  simple  case  of  a  half-wave  linear  rectifier  fed  with  a  single 
sinewave,  and  then  with  the  sum  of  two  sinewaves.  After  this  introduction  we  go  on 
to  the  accelerating  and  compressive  rectifiers  fed  with  the  sum  of  two  sinewaves,  and  fi¬ 
nally  discuss  cascaded  rectifers  and  parallel-cascaded  rectifiers  of  the  same  type  and  of 
mixed  types. 

[1]  HALF-WAVE  LINEAR  RECTIFIER  :  RESPONSE  TO  A  SINGLE  SINUSOID. 

Let  the  input  to  a  half-wave  rectifier  (y  —  ex,  x  >  0;  y  =  0,  x  <  0)  be 
e{t)  =  A  cos(pt  +  0p)  =  A  cos  x,  where  p  =  27r*frequency  of  input  and  0p  —  phase.  Taking 
A  >  0  and  the  constant  of  proportionality  c  =  1,  the  output  is  a  function  f{x),  mwhere 


f[x)  = 


A  cos  x,  cos  x  >  O 
0,  cos  x  <  0. 

We  can  express  f{x)  in  terms  of  a  Fourier  series  in  x,  where 


and 


and 


/(x)  =  a0/2  +  Y2  ancos  nx 


n=l 


1  f  1 

an  =  ~  f(x)cosnxdx ,  /(x)  =  0,  |x|  >  %/2,  n  —  0,1,2,.. 

J-k/2 

_  i 

n  J -ie/2 
A  I"*/2 

=  —  /  [cos(n  +  l)x  ■+■  cos(n  —  l)xl  dx,  n  ^  1 

K  Jo 


A  cos  x  cos  nx  dx 


2A( _ 1\(3+™)/3 

— V  V  ,', -  n  even 


0 


n  odd,  u^l, 


a  i 


=  2/ 
7T  Jo 


ir/2  ^ 

cos 2  xdx~  - 
2 


.  A  A  2A  2 A 

=>  f[x)  = 1 cos  x  H - cos2x - cos  4.x -\ . 

'  7T  2  37T  157T 

(2 j  HALF-WAVE  LINEAR  RECTIFIER  :  RESPONSE  TO  THE  SUM  OF  TWO 
SINUSOIDS 
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If  the  input  voltage  is  given  by 


e(f)  =  Pcos(pt  +  dp)  +  Qcos(qt  +  0q) 
then  we  can  rewrite  this  as 

e(t)  =  P[cos(pt  +  Op)  +  kcos(qt  +  0q)] 

where  k  =  Q/P. 

The  case  k  <  1 

Without  loss  of  generality,  we  can  take  P  >  0  and  the  constant  of  proportionality,  c, 
to  be  1.  First  let  us  consider  k  <  1,  and  set 

y)  -  /  P{cosx  +  keosy),  {cosx +  k  cosy)  >  0 
'  1 0,  (cos x  +  kcosy)  <0 

where  x  =  (pt  +  0p ),  and  y  =  [qt  +  0q). 

f{x,y)  is  a  surface  in  and  above  the  (x,y)-plane,  bounded  by  (cosx  +  kcosy)  =0 
in  the  (x,  y)-plane.  Clearly  adding  2ir  to  x  or  y  leaves  f[x,y)  unaltered,  so  f{x,y)  is  a 
periodic  function  in  x  and  y.  So  if  we  know  /(x,y)  in  the  rectangle  (— 7r,7r)  *  (— 7r,7r)  we 
will  know  all  its  values. 

Since  /(x,  y)  is  bounded  in  the  rectangle  (— tt,  7r)  *  (— 7r,  tt)  and  its  first  derivatives  are 
bounded,  the  double  Fourier  series  in  (x,  y)  of  /(x,y)is  a  valid  expansion  in  this  rectangle 
(Hobson,  1926).  If  the  Fourier  series  of  /(x,y)is  valid  in  the  (x,y)  plane,  then  it  is  valid 
on  the  line  py  —  qx  =  pOq  —  qOp,  found  by  eliminating  t  from  x  =  (pt  +  0P),  y  =  ( qt  +  0q). 

The  boundaries  of  f(x,y) axe  the  curves  given  by  ccsx+  kcosy  =  0.  In  the  shaded 
area,  cos  x  +  kcosy  >  0,  elsewhere  cosx  +  kcosy  <  0,  giving  /(x,y)  =  0.  Since  /(x,y) is 
an  even  function,  its  double  Fourier  expansion  will  be  a  cosine  series  given  by 

f{x>y)  =  r^oo  +  ^2  A±mncos(mx  ±  ny)  +  Aiqcqs  x  +  Aqicos  y 


m=  1  n=  1 


where 


A±mn  =  J  j  f{x,y)cos(mx±ny)dxdy 
i  fw  r 

=  - J  J  f(x,y)\cosmxcosny  sinmxsinny\dxdy. 
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Since  the  region  is  symmetrical  in  both  x  and  y,  A±mn  can  be  found  by  using  one  quarter 
of  the  plane.  Hence 


2  p  fit  rarccoa{—kco3y) 

A±mn  = —r  /  cos  ny  /  (cos  x  +  k  cos  y)cosmx  dx  dy 

*  Jo  Jo 

since  f(x,y )  =  0  when  x  >  arccos(—k  cosy). 

The  calculation  for  A±mn,  when  m  =  2  and  n  =  0,  is  shown  below. 


•'4.20 


tt  i.arccoa{  —  kcoay) 


2P  r  r 

= — 2"  /  /  (cos  x  +  k  cosy)  cos  2x  dx  dy 

n  Jo  Jo 

=  ^~2  Jo  -  k2cos2y)z/2dy 
=  4P  fl(l-k2z2f2  ^ 

3  Wo 

4 P  (  [l(l  ~k2z2\k,  f1  l2  2fl-k2z2\h  ,  1 

=  ^U  d*-]0  k‘*: t-np-)  4 


Using  the  identity 

/  (1  —  k2z2)  \  ^  _  f  1  k2z2  ) 

V  (1  -  z2)  )  l  {(1  —  ib3«2)(l  —  *a)>*  {(1-  k2z2)(l-  z2)}i  / 

and  letting 

- -  dz 

k2z 2)}5 

then  Zo  =  K,  the  complete  elliptic  integral  of  the  first  kind,  and  Z3  can  be  expressed  in 
terms  of  Z3- 2  and  Za-A  by  using  the  recurrence  formula 


-f 


7  _  (s  —  2) (1  +  k2)Z3- 2  —  (s  —  3)Z3_4 
(s  -  1  )k2 


for  s  >  4,  (Bennett,  1933).  From 


Z2  =  (K-E)/k2, 
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where  E  is  the  complete  elliptic  integral  of  the  second  kind,  we  have  that 
A20  =  ^[E  -k2Z2  +  k4Z4\ 

4  p 

=  ^[E  -  {K  -  E)  +  { 2  +  k2)K/ 3  -  2(1  +  k2)E/ 3] 

4  P 

=  Sj[2(2-t2)£-(l-fc2)ir], 

This  gives  the  amplitude  of  the  frequency  ( mx  ±  ny)/2ir  and  the  phase  angle  ( mdp  ± 
nOq).  The  values  of  the  amplitudes  for  m  and  n  =  0, 1,2, 3,4  are  as  follows  : 

4  P 

A00  =  —  [2E  -  (1  -  k2)K\ 

7T 

,  P 

Al0=2 

4  kP 

Aoi  =  — 

4  p 

A20  =  ^[2(2  -  k>)E  -  (1  - 
4  p 

All  =  3^)fcl(1  +  *  )X1 

4  P 

A°2  =  9^F[2(2/c2  ~  ^  +  (2  “  3^)(1  "  *  )*1 
4P 

^4°  =  £^2  [(-38  +  88/t2  -  48 k4)E  +  (23  -  47 fc2  +  24A:4)JFs:] 

4  p 

Mi  =  ^r^;[(8A:4  -  13fc2  +  3)-B  -  (1  -  fc2)(3  -  4fc2)if] 

4P 

Au  =  T^l(k'  ~ " 2)K  ~ 2{k4  -e  +  1)£l 

=  tt~TT3  [{8  —  1342  -f-  3 k')E  -  (8  -  17 42  +  94") JT] 


457T2 /c3 
4P 


2404  =  “  1) ( — 15/c4  +  64A:2  -  48 )K  -  (38 A:4  -  88A:2  +  48)E}. 

The  third  and  higher  odd  order  terms  are  zero,  and 


r  ir/2 

if  =  /  (1  —  k2sin20)~%  d9 

Jo 

=  f  {{l-z2){l-k2z2)]-Uz 

Jo 
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and 


fir/ 1 

E  =  (l  —  k? sin2 0)*  dz 

Jo 


(1  —  k2z2) 2  (1  —  z2)~ 3  dz 

JO 

where  k  <  1. 

The  case  k  >  l. 

We  can  rewrite  f[x,y)  in  the  following  way: 

P(cos  y  +  l  cos  x)/l,  cosy  +  l  cost  >0 
cos  y  -f-  /  cos  x  <  0 

where  /  =  l/k  <  1  and  consequently 

OO  CO 

f{x,y)  =  A'00/2  +  Y]  A'±Tacos(ry  ±  sx)  +  A'10cosy  +•  A'0lcosx 


f[x 


r=l 3=1 


where 


A>  =  *L 

±r3  In2 


r  r 

/  cos  sx 
Jo  Jo 


*arccos(— Icoax) 

cos  sx  I  (cos  y  +  l  cos  x)cos  ry  dy  dx. 


A'^s  is  the  coefficient  of  cos(rx  ±  sx)  which  may  be  written  as  cos(sx  ±  ry).  So  for  a 
given  m  and  n,  say  M  and  N,  we  will  have  to  consider  A±mn>  for  k  <  1  and  A'±NM  for 
k  >  1.  For  example,  let  us  consider  the  coefficient  of  cos  2x. 


^±02  — 


2 P_ 

In 2 
4  P 
9n2l2 
4  Pk3 

9n2 

4P 


fie  farccos{—lco3x) 

/  cos  2x  /  (cos  y  + 1  cos  x)  dy  dx 

Jo  Jo 

[2(2 12  -  1)JS  +  (2  -  3f2)  (1  -  l2)K] 

[2(2/it2  -  l)E(l/k)  +  (2-  3/k2)(l  -  l/k2)K(l/k)] 


~  9lT2k  [2^(2  ~  k2)E(l/k)  +  (2k2  -  3)(k2  -  l)K(l/k)}. 
Therefore  the  function  of  amplitude  g(k)±mn  is  given  by 


g(k) 


±mn 


_  f  A±mn,  k  <  l 

l  A±nm,  k  >  1. 


l±nm 
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4  P 

A'i00  =  -  (k*  -  1)K] 

u  v> 


A'  -kP 

A±1°  -  T 

A'  -P 
±01  ~  2 

4  P 

^±ii  =  ^[(A2  +  i)£-(fc2-i)*:] 

4  p 

A'±™  =  Me  t2*2(2  ~  k2)E  +  ( 2k 2  "  3)^2  ~  l)K) 

^±40  =  22^2 fc3  [(23fc4  ~  47&2  +  24)^  +  (-38*4  +  88A:2  -  48 )E] 

,  4P 

■4±31  =  45^  [<8  - m  +  n')E  -  (k2  -  0  (“2  -  4)JST] 

4  p 

A± 22  =  [(!  -  fc2)(1  -  %k2)K  -  2(1  -  A:2  +  A:4)£] 

4  p 

^±13  =  ^  [(8fc<  ~  13*2  +  3)£  ~  (8*4  ~  1 7k2  +  9)^T] 

4P 

^±04  =  ^2 1  [(X  -  fc2)(-13  +  64A:2  -  48A:4)ir  -  k2{ 38  -  88/t2  + 


48  k*)E). 


The  third  and  higher  odd  order  terms  are  zero,  and  E  and  K  are  functions  of  1  fk  <  1. 
The  elliptical  integrals  were  calculated  using  well-known  algorithms  (King,  1924.  Re<*an 
1985). 

[3J  HALF-WAVE  SQUARE  LAW  RECTIFIER  :  RESPONSE  TO  THE  SUM  OF 
TWO  SINUSOIDS. 

If  the  rectifier  is  of  the  form  y  =  cx2,  x  >  0  and  y  =  0,  x  <  0  and  if  A:  <  1  then,  as  for 
the  half-wave  linear  rectifier,  we  can  consider  the  rectifier’s  output  as  the  function  f{x,y) 
where 


.  f  P2 (cos  X  +  k  cos  y)2 , 

/(*.!,)  =  |  0, 


cos  x  +  k  cos  y  >  0 
cos  x  +  Ar  cos  y  <  0 


where  x  —  (pt  +  9p),  and  y  =  [qt  +  9q).  Again  f(x,y)  is  bounded  in  the  rectangle 
(-7r,7r)  *  (— 7r,?r)  by  cosx  +  kcosy  and  its  Fourier  expansion' will  be  a  cosine  series  given 
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by 

^  03  CO 

f(x,y)  =  -Aoo  +  Yl  A±rnnCOs(mx  ±  ny)  +  Aiqcos  x  +  Apicosy 
m  — 1  n=l 


but  now 


Ai, 


_  2 p2  r 

in  —  o  / 

7T2  Jo 


cos  ny 


parccos(-kcoay) 

/  (cos  x  +  k  cos  y)2 cos  mx  dx  dy 

Jo 


since  f(x,y )  —  0  when  x  >  arccos(—k  cosy).  When  k  >  1,  we  have 

■a  rccosf  —  lcoax) 


±ra  ~  £2^2 


A!  - 


rir  pa 

/  COS  SX 

Jo  Jo 


COS  SX 

2p2^.2  rir  fa.rccoa(—lco3x) 

=  — ^2 —  J  cossx  J  (cos  y  +  l  cos  x)2 cos  ry  dy  dx 


(cos  y  +  l  cos  x)  cos  ry  dy  dx 
"arccoaf—lcoax) 


where  l  =  1/A:. 

[4]  HALF-WAVE  SQUARE-ROOT  RECTIFIER  :  RESPONSE  TO  THE  SUM  OF 
TWO  SINUSOIDS. 

Now  the  rectifier  is  of  the  form  y  =  Cy/x,  x  >  0  and  y  —  0,  x  <  0  and  for  k  <  1  we  will 
have  the  function 


f(x,y) 


P?  (cos  x  +  kcosy)? ,  cosx  +  kcosy>0 
0,  cos  x  +  k  cos  y  <  0 


where  x  =  (pt  +  0P),  and  y  =  (qt  +  9q).  Thus 


2  pi  c*  parccoaf—kcoay) 

■A+mn  =  — 5—  /  cosny  /  (cos  x  +  k  cosy)*  cos  mxdx  dy 

n2  Jo  Jo 

since  f(x,y)  =  0  when  x  >  arccos(—k cosy)  and  for  k  >  1,  we  have 


2Pa  f v  rarccoa(—lco3X)  ^ 

A^rs  =  —£ —  /  cossx  /  (cosy  +  l  cos x)*  cosry  dy  dx 

hn2  Jo  Jo 


parccosf—lcosx) 


2  pkk*  f K  parccoaf—lcosz) 

=  - - —  /  cos  sx  /  (cos  y  +  l  cos  x) 1  cos  ry  dy  dx 

n2  Jo 


Jo 


where  l  =  l/k.  Similarly,  we  can  find  the  response  to  any  half-wave  rectifier  whose  equation 
is  y  =  cxn,  x  >  0;  y  =  0,  x  <0,  where  n  is  any  real  number. 
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[5]  TWO  CASCADED  LINEAR  HALF-WAVE  RECTIFIERS,  A.C.  COUPLED. 

If  two  rectifiers  are  D.C.  coupled,  the  output  will  be  the  same  as  a  single  linear  half¬ 
wave  rectifier.  Indeed,  if  two  half-wave  rectifiers  are  D.C.  coupled  and  the  first  of  the  series 
is  a  linear  rectifier,  the  final  output  will  be  the  same  as  that  of  the  second  rectifier  alone. 
After  the  two  sinusoids  pass  through  the  first  rectifier,  their  function  is  given  by 


(  P(cos  x  +  k  cos  y),  (cos  x  +  k  cosy)  >  0 
f(x,  y)  =  < 

[0,  (cos  x  +  k  cos  y)  <  0 

where  x  =  (pt  +  0p),  and  y  —  (qt  +  0q).  This  has  a  D.C.-level  given  by  Aoo/2,  the  constant 
term  in  the  double  Fourier  series  expansion  of  f(x,y).  If  our  two  successive  rectifers  are 
linked  by  A.C.  coupling,  this  D.C.-level  must  be  removed  and  so  the  function  entering  the 
second  rectifier  is  given  by 

F(x,y)  =  f(x,y)  -  A00/2 


where 


arcco3{—kcosy) 


( cos  x  +  k  cos  y)  dx  dy. 


After  passing  through  the  second  rectifier,  the  output  is  given  by 


<f>(x>y) 


F(x,y),  F(x,  y)  >  0 
0,  F(x,  y)  <  0. 


This  can  be  represented  by  a  double  Fourier  series  where  the  coefficients  A±mn  are  given 
by 

2 p  r*  r* 

A±mn  =  ~ r  /  cosny  /  <f>(x,y)cosmxdxdy. 

*  Jo  Jo 

[6]  CASCADED  COMPRESSIVE  RECTIFIERS 

Results  have  been  obtained  for  two  square  root  ( y  =  cx?,x  >  0;y  =  0,z  <  0) 
rectifiers  in  series  in  particular,  for  any  combination  of  two  compressive  rectifiers  (y  = 
cxN  ,x  >  Q;y  =  0,  x  <  Oand  <  N  <  1)  and  for  three  square  root  rectifiers  in  series. 

(7 j  TWO  PARALLEL  LINEAR  RECTIFIERS  WHOSE  SUMMED  OUTPUTS 
PASS  THROUGH  A  THIRD  LINEAR  RECTIFIER  :  THE  DICHOPTIC  OR  DI- 
CHOTIC  CASE 
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In  this  situation  one  only  frequency  (FT)  passes  through  rectifier  no.  1  and  only  one 
frequency  (.F2)  passes  through  rectifier  no. 2  in  parallel  with  the  first  rectifier.  Then  the 
output  from  both  rectifiers  combine  to  form  the  input  of  the  third  rectifier. 

The  output  of  the  first  rectifier  is  f(x)  where 


/(*)  = 


P  COS  X,  COS  X  >  O 
0,  cos  x  <  0 


with  a  D.C.-level  of  P fit.  The  output  of  the  second  rectifier  is  g{y )  where 


g{y)  = 


Pk  cos  y,  cos  y  >  O 
0,  cos  y  <  0 


whose  D.C.-lcvel  is  Pk/ n.  To  adjust  for  the  D.C.-level,  the  input  to  the  third  rectifier  will 
be  the  function 

h(x,y)  =  f(x)  -  P/ 7r  +  g(y)  -  Pk/ir. 

The  output  from  the  third  rectifier  is  given  by 


-)={o! 


h{x,y),  h(x,y)>  0 
0,  h(x,y)  <  0. 


Hence  the  coefficients  of  the  double  Fourier  series  can  be  found  for 


2 p  [K  fw 

A±mn  =  —rr  /  cosny  /  H(x,y)cosmxdxdy 

Jo  Jo 

Results  have  been  obtained  for  the  case  that  all  three  rectifiers  have  a  linear  characteristic 
and  coupling  is  A.C.  rather  than  D.C.  Other  cases  such  as  mixed  rectifiers  (e.g.  where 
nos.l  and  2  are  cube  root  rectifiers  and  no.3  is  a  square  law  rectifier)  are  amenable  to  the 
same  general  mathematical  treatment. 

(8)  HALF-WAVE  RECTIFIER  COMBINING  ACCELERATING  AND 
COMPRESSIVE  SEGMENTS 

For  this  rectifier,  the  curve  equation  is  given  by 


{0,  x  <  c 

d(x  —  c)‘l,  c  <  x  <  5c 

(x  —  c) 1//1<3  —  g,  5 c  <  x 
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where  d  =  1/64(4c)tI  and  g  =  || (4c)  i«  and  c  is  chosen  suitably.  Consequently 

!0,  cos  2  +  kcos  y  <  c 

P4d(cos  x  +  kcosy  -  |>)4,  c  <  cos  x  +  kcosy  <  5c 

P To  (cos  x  +  k  cos  y  —  —  g,  5c  <  cos  x  +  k  cos  y 

where  x  =  (pt  +  0p),  and  y  =  (qt  +  0q).  So 


A±mn  =  -t;  [  cosny  [  f  (x,  y)  cos  mx  dx  dy 
ft1  Jo  J o 
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zoom-FF 

According  to  the  Heisenberg-Gabor  uncertainty  principle  the  limiting  frequency  resolution 
of  a  spectrum,  AF  Hz,  is  given  by 

AF  =  1/AT 

where  AT  is  the  recording  duration.  Thus,  for  example,  a  recording  of  duration  500  sec  could,  in 
principle,  be  analyzed  at  a  resolution  of  0.002  Hz  so  that,  if  the  bandwidth  were  DC-100  Hz,  the 
spectrum  would  contain  100  x  500  =  50,000  lines.  In  practice,  however,  the  FFT  usually 
provides  many  fewer  lines,  typically  several  hundred  over  a  DC-100  Hz  bandwidth.  We  have 
developed  a  nondestructive  form  of  zoom  FFT  that  allows  high  zoom  ratios  (typically  32-64)  over 
a  wide  bandwidth  so  that  we  routinely  obtain  25,000  or  50,000  lines  over  DC-100  Hz. 

The  method  is  to  digitize  a  time  series  of  duration  AT  by  means  of  a  Bruel  and  Kjaer 
spectral  analyzer.  The  digitized  time  series  is  recorded  on  floppy  disk  in  an  Hewlett-Packard 
model  9000  computer  that  controls  the  analyzer.  If,  for  example,  the  bandwidth  is  DC-100  Hz, 
the  sampling  rate  will  be  250  Hz.  We  routinely  digitize  a  320-sec  duration  of  the  time  series. 
Next,  the  digitized  data  are  replayed  at  much  increased  rate  (25  kHz  rather  than  250  Hz),  filtered 
and,  for  example,  the  DC-3.0  Hz  section  submitted  to  FFT,  giving  800  lines  within  DC-3.0  Hz. 
This  destroys  the  time  series  in  the  analyzer.  Now  the  time  series  is  replayed  again  at  25  kHz, 
heterodyned  to  shift  the  3.0-6.0  Hz  segment  to  DC-3.0  Hz,  filtered,  resampled,  submitted  to  FFT 
and  shifted  back  to  3.0-6.0  Hz.  This  gives  us  800  lines  within  3.0-6.0  Hz.  The  process  is 
repeated  to  give  800  lines  in  each  3.0  Hz  segment  between  DC  and  100  Hz. 

The  value  of  this  method  in  electrophysiology  is  not  self-evident.  The  value  is  based  on 
our  fortunate  discovery,  illustrated  in  Figure  15,  that  the  discrete  frequency  components  of  the 
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steady-state  evoked  potential  are  of  ultra-narrow  bandwidth,  and  can  be  less  than  0.002  Hz. 
Consequently,  the  noise  is  spread  through  50,000  bins  while  signal  components  are  concentrated 
into  one  or  two  bins.  This  gives  (a)  high  signal-to-noise  ratios,  and  (b)  excellent  separation  of 
signal  components.  The  procedure  has  been  published  in  a  book,  sponsored  in  part  by  AFOSR(^) 
and  is  also  in  press  in  a  journal  article.^) 


C\J 

0 


1 .0 


.8 


Z-  • G 


O 

> 

w 

L. 

0) 

3 

o 

CL 


.4 


.2 


0.0 


L 


31. B0  31.85 

F  requency  ( Hz ) 


Figure  15.  Two-Thousandth  of  a  Hertz  Bandwidth  At  a  frequency  of  0.0019  Hz 
approximately  two  thirds  of  the  EP  signal  power  is  contained  within  a  single  bin. 
Therefore  the  EP  signal's  bandwidth  is  no  more  than  0.0019  Hz.  Note  that  even  if 
the  EP  signal's  bandwidth  were  considerably  less  than  0.0019  Hz,  more  than  one 
bin  of  the  analyzer  would  contain  power  when  the  analyzer's  resolution  is  set  at 
0.0019  Hz,  because  he  sensitivity  of  adjacent  bins  overlaps.  The  stimulus  wa  a 
vertical  sine  wave  grating  that  was  counterphase-modulated  at  nominally  F  =  8 
Hz.  The  4F  component  is  shown.  Zoom  FFT  was  applied  to  the  entire  length  of 
EEG  samples  of  520-sec  duration.  To  allow  some  frequency-domain  averaging 
( with  overlap ),  recording  durtion  was  640  sec. 
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(11)  Use  of  the  two-sinewave  method  to  measure  orientation  tuning  in  human  cortical  neurons 

A  vertical  sinewave  grating  of  spatial  frequency  5  c/deg  was  generated  on  a  Joyce  CRT  and 
counterphase-modulated  at  frequency  FI  (nominally  8  Hz).  A  second  grating  of  spatial  frequency 
5.5  c/deg  and  variable  orientation  was  generated  on  a  second  Joyce  CRT  and  counterphase- 
modulated  at  frequency  F2  (nominally  7  Hz).  The  two  gratings  were  optically  superimposed. 
Field  size  was  10  deg,  contrast  was  40%  for  each  grating  and  mean  luminance  was  250  cd/m2. 
Calibration  with  a  linear  photocell  showed  that  each  CRT  was  quite  linear:  second  harmonic 
distortion  was  below  0.1%  of  the  fundamental  component's  power.  Cross-modulation  terms  were 
essentially  zero  because  different  CRTs  driven  by  different  electronics  generated  the  FI  Hz  andF2 
Hz  gratings.  Photocell  calibration  showed  cross-modulation  components  to  be  less  than  0.01%  of 
the  fundamental  components'  power. 

Human  steady-state  evoked  potentials  were  recorded  between  electrodes  placed  on  the  inion 
and  midway  between  the  inion  and  the  vertex  along  the  midline.  Responses  were  analyzed  in  the 
frequency  domain  by  a  Bruel  and  Kjaer  analyzer  (model  2032)  modified  to  carry  out  zoom-FFT 
nondestructively  at  high  zoom  factors  over  a  wide  bandwidth. (29)  Resolution  was  0.0078  Hz  over 
a  DC-100  Hz  bandwidth  for  a  320-sec  recording  period,  i.e.  12,800  frequency  bins  were  available 
with  frequency-domain  averaging  also. 

The  dashed  line  in  Figure  16  plots  the  amplitude  of  a  (2F1  +  2F2)  cross-modulation  res¬ 
ponse  term  as  a  function  of  the  orientation  difference  between  the  gratings.  This  cross-modulation 
term  necessarily  indicates  a  nonlinear  interaction  between  visual  responses  to  the  fixed  vertical 
grating  and  the  variable-orientation  grating,  and  has  previously  been  shown  to  be  substantially 
independent  of  spatial  phase.(31.32)  Figure  16  shows  that  the  nonlinear  interaction  was  large  when 
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the  gratings  were  parallel  and  fell  to  a  minimum  when  their  orientations  differed  by  about  30  deg. 
The  half-height  full  bandwidth  of  the  curve  is  about  12  deg.  The  frequency-doubled  2F1  Hz 
response  produced  by  the  fixed  vertical  grating  was  suppressed  when  the  two  gratings  were 
parallel,  but  the  second  grating  had  comparatively  little  effect  when  grating  orientations  differed  by 
about  30  deg.  Similar  results  were  obtained  for  a  second  and  third  subject. 


Figure  16.  Nonlinear  interactions  between 
responses  to  two  gratings  as  a  function  of 
orientation  difference.A  vertical  grating  was 
counterphase-modulated  at  FI  Hz  and  a 
superimposed  variable-orientation  grating 
was  modulated  at  F2  Hz.Solid  symbols  plot 
the  amplitude  of  the  nonlinear  cross-modu¬ 
lation  (2 FI  +  2F2)Hz  term  in  the  evoked 
potential  versus  the  variable  grating's 
orientation.  Open  symbols  plot  the 
frequency-doubled  2F1  Hz  term.  Results 
are  shown  for  two  subjects. 
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The  observations  reported  above  can  be  understood  if  the  (2F1  +  2F2)  term  is  generated  by 
cortical  neurons  tuned  to  a  narrow  range  of  orientations  (such  as  those  described  by  De  Valois  et 
al.W  When  the  grating  orientations  differ  by  more  than  about  30  deg,  most  of  these  neurons 
cannot  encompass  both  gratings  within  their  orientation  bandwidths,  and  will  therefore  fail  to 
generate  cross-modulation  terms. 

However,  when  we  placed  the  two  gratings  at  right  angles  (the  fixed  grating  remaining 
vertical),  the  nonlinear  cross-modulation  term  rose  to  a  second  maximum.  For  subject  B  this  (2F1 
+  2F2)  term  was  as  large  for  near-orthogonal  gratings  as  for  parallel  gratings,  and  only  a  little  less 
for  subject  A.  The  interaction  term  was  largest  at  exactly  90  deg  orientation  difference  for  subject 
A  but,  curiously,  peaked  sharply  just  5  deg  from  90  deg  for  subject  B. 

This  finding  that  there  is  a  strong  nonlinear  interaction  between  responses  to  vertical  and 
near-horizontal  gratings  can  be  understood  if  we  assume  that  cortical  neurons  tuned  to  a  narrow 
range  of  orientations  around  the  vertical  interact  nonlinearly  with  cortical  neurons  tuned  to  a  narrow 
range  of  orientations  around  the  horizontal.  It  may  be  relevant  that  conical  neurons  tuned  to 
different  orientations  can  inhibit  each  other  when  excited  simultaneously/^4-55) 

If  our  findings  can  be  generalized  to  other  kinds  of  two-dimensional  pattern,  this  would 
imply  that  human  VEPs  to  patterns  modulated  in  two  dimensions  cannot  entirely  be  explained  in 
terms  of  VEPs  to  gratings.  In  particular,  the  findings  reported  here  could  not  result  from  the 
stimulation  of  independent,  linear,  orientation- selective  mechanisms. 

(12)  Nonlinear  interactions  between  visual  responses  to  two  gratings  of  different  spatial 
frequencies 
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We  have  measured,  interactions  between  responses  to  two  parallel  or  orthogonal  gratings  of 
different  spatial  frequencies.  Subjects  viewed  a  fixed  1.0  c/deg  vertical  grating  counterphase 
modulated  at  Fi  =  8  Hz  superimposed  on  a  grating  of  variable  spatial  frequency  and  orientation, 
modulated  at  F2  =  7  Hz.  We  analysed  brain  responses  using  high-resolution  zoom-FFT,  at  a 
frequency  resolution  of  0.008  Hz  (12,500  lines  DC-100  Hz).  Nonlinear  processing  produces 
multiple  discrete  components  of  frequency  (nFi  ±  ni?2)  where  n,  m  are  zero  or  integral.  We 
recorded  up  to  20  such  discrete  components,  each  with  a  bandwidth  as  narrow  as  0.002  Hz.  With 
parallel  vertical  gratings,  the  2Fi  Hz  response  of  the  1  cl deg  grating  was  suppressed  nonlinearly 
by  a  second  grating  of  higher  spatial  frequency  over  a  broad  1.5-13  c/deg  range,  but  the 
suppression  did  not  occur  in  the  opposite  direction  from  low  to  high  spatial  frequency. 
Interactions  between  orthogonal  gratings  were  qualitatively  different.  The  2Fi  Hz  response  of  the 
vertical  1.0  c/deg  grating  and  the  2F2  Hz  response  of  a  1.5-4  c/deg  horizontal  grating  were  both 
attenuated,  suppression  being  exerted  both  ways  along  the  spatial  frequency  axis.  And  with  a 
horizontal  grating  of  4-13  c/deg  there  was  comparatively  little  suppression  of  2Fi  and  2F2 
components. 

This  work  was  reported  at  ARVO(36)  and  a  manuscript  is  in  preparation. 

(13)  Animal  experiments:  Nonlinear  analysis  of  single-cell  firing  and  local  slow  potentials  in 
anesthetized  cat. 

This  is  a  joint  study  (M.P.  Regan,  K.  Grasse  and  D.  Regan)  and  is  in  progress.  The  aim  is 
to  apply  the  mathematical  and  zoom-FFT  approach  described  in  Section  2b  (4)  above  to  single-cell 
recording  of  cortical  cells  so  as  (A)  to  develop  a  means  of  modelling  the  multi-neuron  sequence  of 
processing  between  retina  and  cortex,  (B)  to  develop  a  new  means  of  studying  nonlinearities  of 
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single  cells,  and  (C)  to  find  whether  the  very  considerable  stability  of  neural  signals  implied  by  the 
narrow  bandwidth  of  the  human  evoked  potential  signals  (Figure  13)  is  evident  at  single  cell  level 
or  whether  it  is  a  multi-neuron  property. 

In  order  to  carry  out  this  study  I  have  developed  an  electronic  device  for  converting  the 
instantaneous  firing  frequency  of  a  single  neuron  into  an  analog  signal  that  can  be  subjected  to 
nondestructive  zoom-FFT  analysis  as  described  in  Section  2b  (4)  above.  Visual  stimulation 
equipment  has  been  installed  in  Dr.  Grasse's  animal  laboratory  and  connected  to  the  zoom-FFT 
analyser  in  my  laboratory.  We  have  recorded  from  three  cats  and  shown  that  nonlinear  interactions 
between  responses  to  two  gratings  can  be  analysed  in  cortical  slow  waves.  We  have  recorded 
single-cell  firing  in  striate  cortex,  and  established  the  appropriate  stimulus  conditions. 

(14)  Installation  of  the  BTi  7-channel  Neuromagnetometer  and  magnetically  shielded  room 

Installation  of  the  magnetically  shielded  room  started  on  September  12,  1988  in  a  room  set 
aside  for  the  purpose  in  the  Farquarson  Building  at  York  University.  Installation  was  completed 
on  schedule.  Installation  of  the  neuromagnetometer  was  started  on  October  5,  1988.  The  system 
dewar  was  cooled  to  liquid  helium  temperature  during  the  week  of  November  7,  and  has  been 
maintained  at  liquid  helium  temperatures  since  then.  Several  of  the  students,  technicians  and 
faculty  at  York  have  been  trained  to  transfer  liquid  helium  from  a  100-litre  reservoir  to  the  system 
dewar.  This  must  be  done  three  times  a  week.  The  total  expenditure  of  liquid  helium  is  stabilizing 
at  about  100  litres  per  two  weeks.  BTi  representatives  continued  to  install  and  check  out  the 
magnetometer  up  to  November  11,  and  throughout  the  week  of  November  14-18,  five  of  us  were 
instructed  by  BTi  representatives  on  the  use  of  the  computer  system  and  recording  procedures. 
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In  early  1989  BTi  installed  a  sensory  position  indicator  (SPI)  on  the  magnetometer. 
Subsequent  to  installation  one  of  the  recording  channels  developed  an  intermittent  fault.  After 
considerable  investigation  this  was  shown  in  November  1989  to  be  caused  by  hardware  damage, 
probably  caused  during  SPI  installation. 

(15)  Analysis  of  human  magnetic  brain  responses  by  nondestructive  zoom-FFT 

Figure  17  shows  evoked  magnetic  brain  responses  analysed  by  zoom-FFT  at  a  resolution 
of  0.008  Hz.  In  Figure  17A  the  subject's  contralateral  ear  was  stimulated  with  a  1000  Hz  tone  that 
was  amplitude  modulated  at  Fi  =  20.351  Hz.  The  spike  was  at  2Fi  =  40.703  Hz.  In  panel  B  the 
1000  Hz  tone  was  modulated  at  F2  =  20.742  Hz.  This  spike  was  at  2F2  =  41.484  Hz.  In  panel 
C,  the  ear  was  stimulated  with  both  AM  tones  simultaneously.  The  spikes  in  A  and  B  were 
suppressed,  and  a  new  spike  at  (Fi  +  F2)  =  41.093  Hz  emerged. 

The  very  narrow  bandwidth  of  the  spikes  indicates  that  the  steady-state  magnetic  response 
is  very  stable  in  amplitude  and  phase.  This  was  observed  at  all  recording  sites  (only  one  is 
shown).  The  nonlinear  phenomenon  shown  in  Figure  17  was  replicated  when  we  applied  the 
stimulus  to  a  hardware  model  of  the  ear’s  hair  cell  transducer  function. 

One  of  use  (M.P.  Regan)  is  developing  a  mathematical  approach  to  modelling  responses 
such  as  those  shown  in  Figure  17  so  as  to  allow  neuromagnetic  data  to  be  used  to  test  multi-neuron 
models  of  the  auditory  pathway. 

This  work  is  in  press(37>18)  and  a  review  article  has  been  published^39). 

( 1 6)  Joint  work  with  Professor  L.  Kaufman 

A  joint  experiment  with  Profesor  L.  Kaufman  on  the  effects  of  attention  on  visual 
responses  has  been  planned,  and  we  have  constructed  the  required  apparatus. 
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Figure  17.  A  small  section  of  the  MEG  spectrum  recorded  during  auditory 
stimulation.  Analysis  was  by  Zoom-FFT  at  a  resolution  of  0.008Hz. 
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(17)  Book:  "Human  Brain  Electrophysiology:  Evoked  potentials  and  evoked  magnetic  fields  in 
science  and  medicine"  by  D.  Regan 

Published  by  Elsevier  1989.  This  is  a  single-author  book  whose  writing  was  sponsored  in 
pan  by  AFOSR.  820  pp,  372  figures. 

This  book  attempts  to  link  (1)  our  knowledge  of  evoked  electrical  and  magnetic  responses 
of  the  human  brain  to  (2)  sensory  perception  and  cognition  and  (3)  the  properties  of  single 
neurons  in  primate  brain.  It  covers  vision,  hearing,  somatosensation  and  cognition.  There  are 
three  parts:  technical  and  mathematical  aspects  of  recording  techniques,  basic  research,  and  clinical 
applications.  To  date  three  journals  have  reviewed  the  book,  and  these  reviews  are  quite 
favourable.  Copies  are  enclosed. 

(18)  Editor  of  two  books:  "Binocular  Vision"  and  "Spatial  Form  Vision" 

Macmillan  is  producing  a  series  of  about  14  volumes  under  the  general  title  "Vision  and 
Visual  Abnormalities."  I  was  invited  to  edit  two  of  these  books.  My  aim  was  to  choose  authors 
who  had  at  least  played  an  important  innovative  role  in  the  development  of  their  topic  over  the  last 
10-20  years  and,  preferentially,  initiated  major  advances  in  their  topic.  In  this  way  I  hoped  that 
authors  would  produce  unique  insights  into  how  modem  understanding  of  the  topics  really  did 
emerge  so  as  to  provide  students  with  a  first-hand  understanding  of  creative  science  that  is  often 
lacking  in  second-hand  accounts.  The  authors  were  asked  to  review  their  topics  at  the  level  of  a 
senior  researcher  while  making  the  chapter  accessible  to  graduate  students.  The  teaching  aspect 
was  emphasized. 
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I  was  fortunate  that  almost  all  of  my  first  choice  authors  agreed  to  contribute,  and  only  very 
few  topics  had  to  be  omitted.  Both  books  are  in  press  for  a  1990  publication  date.  I  am  confident 
that  the  books  will  be  of  considerable  use  to  the  psychophysics,  human  factors  and  single-unit 
research  communities. 


MACMILLAN  VOL.  10  "BINOCULAR  VISION" 

H.  Collewijn,  R.M.  Steinman,  C.J.  Erkelens  and  D.  Regan,  "Binocular  fusion,  tereopsis  and 
Stereoacuity  with  a  Moving  Head". 

J.M.  Foley,  "Binocular  spare  perception" 

R.  Fox,  "Binocular  rivalry" 

R.  Held,  "Development  of  binocularity  and  stereopsis" 

A.E.  Kertesz,  "Cyclofusion" 

H.  Ono,  "Binocular  visual  directions  of  an  object  when  seen  as  single  or  double" 

G.  Poggio,  "Physiological  basis  of  stereoscopic  vision" 

D.  Regan,  "Depth  from  motion  and  motion  in  depth" 

C.  Schor,  "Binocular  sensory  disorders" 

C.W.  Tyler,  "The  horopter" 

C.W.  Tyler,  "Cyclopean  vision" 
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MACMILLAN  II  "SPATIAL  VISION” 

J. R.  Bergen,  "Theories  of  visual  texture  perception" 

I.  Bodis-Wollner  and  D.  Regan,  "Spatio-temporal  vision  in  Parkinson's  disease  and  MPTP  treated 
monkeys:  the  role  of  dopamine" 

D.  Levy,  "Spatial  vision  in  amblyopia" 

G.  Mohn  and  J.  Van  Hof-van  Duin,  "Development  of  spatial  vision"  ■ 

M.  Morgan,  "Hyperacuity" 

G. T.  Plant,  "Temporal  properties  of  normal  and  abnormal  spatial  vision" 

D.  Regan,  "A  brief  review  of  some  of  the  stimuli  and  analysis  methods  used  in  spatiotemporal 
vision  research" 

D.  Regan,  "Spatial  vision  for  objects  defined  by  colour  contrast,  binocular  disparity  and  motion 
parallax" 

D.  Regan,  "Spatial  vision  in  multiple  sclerosis" 

K.  Ruddock,  "Spatial  vision  after  cortical  lesions" 

R.  Shapley,  "Contrast  sensitivity:  neural  mechanisms" 

H. R.  Wilson,  "Psychophysical  models  of  spatial  vision  and  hyperacuity" 
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